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Abstract 

We find that in order to completely describe the spin transport, apart from spin current (or linear 
spin current), one has to introduce the angular spin current. The two spin currents respectively 
describe the translational and rotational motion (precession)of a spin. The definitions of these spin 
current densities are given and their physical properties are discussed. Both spin current densities 
appear naturally in the spin continuity equation. Moreover we predict that the angular spin current 
can also induce an electric field E, and in particular E scales as 1/r 2 at large distance r, whereas 
the E field generated from the linear spin current goes as 1/r 3 . 

PACS numbers: 72.25.-b, 85.75.-d, 73.23.-b 
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I. INTRODUCTION 



Recently, a new sub-discipline of condensed matter physics, spintronics, is emerging 
rapidly and generating great interests.— >2 The spin current, the most important physical 
quantity in spintronics, has been extensively studied. Many interesting and fundamental 
phenomena, such as the spin Hall effec t 3l4lSl6 and the spin precession^ in systems with 
spin-orbit coupling, have been discovered and are under further study. 

As for the charge current, the definition of the local charge current density j e (r,t) = 
Re[^(r,t)evty(r,t)] and its continuity equation -^p e (r,t) + V • j e (r,t) = is well-known 
in physics. Here ^f(r,t) is the electronic wave function, v = r is the velocity operator, and 
p e (r,t) = eW^ is the charge density. This continuity equation is the consequence of charge 
invariance, i.e. when an electron moves from one place to another, its charge remains the 
same. However, in the spin transport, there are still a lot of debates over what is the correct 
definition for spin current.— The problem stems from that the spin s*is no invariant quantity 
in the spin transport, so that the conventional defining of the spin current < vs > is no 
conservative. Recently, some studies have begun investigation in this directio n 1 ^ 11 ' 12 , e.g. a 
semi-classical description of the spin continuity equation has been propose d 1 ''' 111 , as well as 
studies introducing a conserved spin current under special circumstances.— 

In this paper, we study the definition of local spin current density. We find that due to 
the spin is vector and it has the translational and rotational motion, one has to use two 
quantities, the linear spin current and the angular spin current, to completely describe the 
spin transport. Here the linear spin current describe the translational motion of a spin, and 
the angular spin current is for the rotational motion. The conventional linear spin current 
has been extensively studied. However, the physical meaning of the angular spin current 
is given for the first time. The definition of two spin current densities are given and they 
appear naturally in the quantum spin continuity equation. Moreover, we predict that the 
angular spin current can generate an electric field similar as with the linear spin current, 
and thus contains physical consequences. 

The paper is organized as follows. In Section II, we first discuss the flow of a classical 
vector. The flow of a quantum spin is investigated in Section III. In Section IV, we study 
the problem of electric fields induced by spin currents. Finally, a brief summary is given in 
Section V. 
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II. THE FLOW OF A CLASSICAL VECTOR 



Before studying the spin current in a quantum system, we first consider the classical 
case. Consider a classical particle having a vector m (e.g. the classical magnetic moment, 
etc.) with its magnitude \m\ fixed under the particle motion. To completely describe this 
vector flow (see Fig.lc), besides the local vector density M(r, t) = p(r,t)m(r,t) one needs 
two quantities: the linear velocity v(r,t) and the angular velocity Cu(r, t). Here p(r,t) is the 
particle density, and v and u describe the translational and rotational motions, respectively 
(see Fig. la and b). In contrast with the flow of a scalar quantity in which one only needs 
one quantity, namely the local velocity v(r,t), to describe the translational motion, it is 
essential to use two quantities, v(r,t) and uj(r,t), for the vector flow. We emphasize that 
it is impossible to use one vector to describe both translational and rotational motions 
altogether. 

Since \rh\ is a constant, the change of the local vector M(r, t) in the volume element 
AV = AxAyAz with the time from t to t + dt can be obtained (see Fig.lc): 



d(M(r,t)Av) = J2 



AV - AV - - 

-—vJr, t)dtM(r, t) - -t-vAt + Ai, t)dtM(v + Ai, t) 
Ai Ai 



i=x,y,z 

+df(r, t) x M(r, t) AVdt (1) 

The first and the second terms on the right describe the classical vector flowing in or out the 
volume element AV and its rotational motion respectively, and both can cause a change in 
the local vector density. When AV goes to zero, we have the vector continuity equation^ 

-rM(r, t) = -V • v(r, t)M(r, t) + cJ(r, t) x M(r, t), (2) 

(JjL 

where vM is a tensor, and its element (vM)ij = ViMj. Notice that this continuity equation is 
from the kinematics and the invariance of \m\, in particular it is independent of the dynamic 
laws.— This is completely same with the continuity equation with a scalar quantity (e.g. 
charge). Introducing j s (r, t) = v(r,t)M(r,t) and j w (r, t) = uj(r,t) x M(r, t), then Eq.(2) 
reduces to: 

j t M{v,t) = -V. j s (r,t) +j u (r,t). (3) 

Here j s = vM is from the translational motion of the classical vector m, and j u = Co x M 
describes its rotational motion.- 13 Since v and uj are called as the linear velocity and the 
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angular velocity respectively, it is natural to name j s and j u as the linear and the angular 
current densities. 

In order to describe a scalar (e.g. charge) flow, one local current j e (r,t) is sufficient. 
Why is it required to introduce two quantities instead of one to describe a vector flow? 
The reason is that the scalar quantity only has the translational motion, but the vector 
quantity has two kinds of motion, the translational and the rotational. So one has to use 
two quantities, the velocity v and the angular velocity uj, to describe the motion of a single 

— * — * — * 

vector. Correspondingly, two quantities j s = vM and j w = uj x M are necessary to describe 
the vector flow. 

In the steady state case, the scalar (e.g. charge) continuity equation reduces into V»j e = 
0, so the scalar current j e is a conserved quantity. But for a vector flow, the linear vector 
current j s is not conserved since V • j s = j w . Whether it is possible to have a conserved 
vector current through redefinition? Of course, this redefined vector current should have 
a clear physical meaning and is measurable. In our opinion, this is almost impossible in 
the 3-dimensional space. The reasons are: (i) One can not use a 3-dimensional vector to 
combine both v and uj. Therefore one can also not use a 3-dimensional tensor to combine 
both j s = vM and j u = uj x M. (ii) Consider an example, as shown in Fig. 3a, a one- 
dimensional classical vector flowing along the x-axis. When x < 0, the vector's direction is 
in +x axis. At < x < L, the vector rotates in accompany with its translational motion. 
When x > L, its direction is along +y axis. Since for x < and x > L the vector has 
no rotational motion, the definition of the vector current is unambiguous, and the non-zero 
element is j xx for x < 0, and j xy for x > L. Therefore, the vector current is obviously 
different for x < and x > L, and the vector current is non-conservative. 

In our opinion, j s = vM and j u = uj x M already have clear physical meanings. They 
also completely and sufficiently describe a vector flow. One may not need need to enforce 
a conserved current. In particular, as shown in the example of Fig. 3a, sometimes it is 
impossible to introduce a conserved current. 

III. THE FLOW OF A QUANTUM SPIN 

Now we study the electronic spin s in the quantum case. Consider an arbitrary wave func- 
tion ^(r, t). The local spin density s at the position r and time t is: s(r, t) = ^(r, f)sM/(r, t), 
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where s — ^a with a being the Pauli matrices. The time- derivative of s(r,t) is: 



d s h 

—siv.t) = - 
dt v ' ; 2 



+ ^e?^-^ }> . (4) 



From the Schrodinger equation, we have 4r^(r,t) = 4i7^(r, t) and 4^(r,t) = 
3Mff\I/(r, £)]'. Notice here the transposition in the symbol f only acts on the spin indexes. 
By using the above two equations, the Eq.(4) changes into: 

(d/dt)s(r, t) = [tfffHV - (H^a^]/2i. (5) 

In the derivation below, we use the following Hamiltonian: 

-2 

H= ^- + V(r) + a*B + ^z»(axp). (6) 

2m n 

Note that our results are independent of this specific choice of the Hamiltonian.— In Eq.(6) 
the 1st and 2nd terms are the kinetic energy and potential energy. The 3rd term is the Zee- 
man energy due to a magnetic field, and the last term is the Rashba spin-orbit couplin g 16 i 17 , 
which has been extensively studied recently^S. Next we substitute the Hamiltonian Eq. (6) 
into Eq.(5), and the Eq.(5) reduces to: 

d - i — > r-r+rP —*\ l — * t "l 

— s = — V • Re{W\— + -(z x a)]aV} 
dt 2 m n 

+ Re{&[B + ^pxz}x (7) 
Introducing a tensor j s (r, t) and a vector j^(r, t): 

j.(r,f) = Re{¥[^ + ^(zxt)]h>} (8) 

l(r,t) = Re{¥^[B + ^{px z)} x JM/}, (9) 



then Eq.(7) reduces to: 

d 



dt s{r,t) = -V«j 8 (r,t)+j w (r,t), (10) 



or it can also be rewritten in the integral form: 

d_ 
dt 




sdV = -fdS*} s + JuJ dV. (11) 
v Js J J Jv 



Due to the fact that = |r = £ + \{z x a) and ±ff = ±[a, H] = \[B + |p x z) x a, Eqs.(8) 
and (9) become: 

j,(r,t) = Re{¥{r,t)dkf{r,t)} (12) 
J w (r,t) = Re{¥(ds/dt)V} = Re{&d x kf}, (13) 



where to = |[5 + ^(p x z)\ is the angular velocity operator.— 

Eq.(lO) is the spin continuity equation, which is same with the classic vector continuity 
equation (3) although the derivation process is very different. In some previous works, this 
equation has also been obtained in the semiclassical case . 1 ^ 11 Here we emphasize that this 
spin continuity equation (10) is the consequence of invariance of the spin magnitude \s\, 
i.e. when an electron makes a motion, either translation or rotation, its spin magnitude 
\s\ = | remains a constant. And the Eq.(10) should be independent with the force (i.e. 
the potential) and the torque, as well the the dynamic law. The two quantities j s (r,t) and 
j w (r, t) in Eq.(10), which are defined in Eqs.(12,13) respectively, describe the translational 
and rotational motion (precession) of a spin at the location r and the time t. They will 
be named the linear and the angular spin current densities accordingly 18 , similar as v and 
uj are called the linear and the angular velocities. In fact, the linear spin current j s (r,t) is 
identical with the conventional spin current investigated in recent studies^. However, we 
give the physical meaning of j w (r, t) for the first time. 

Next, we discuss certain properties of j s (r,t) and j u (r,t). Notice that j w (r,t) which 
describe the rotational motion (precession) of the spin plays a parallel role in comparison 
with the conventional linear spin current j s (r, t) for the spin transport. (1) In similar with the 
classical case, it is necessary to introduce the two quantities j s (r, t) and jaj(r, t) to completely 
describe the motion of a quantum spin. (2) The linear spin current is a tensor. Its element, 
e -g- js,xyi represents an electron moving along the x direction with its spin in the y direction 
(see Fig. 2a). The angular spin current j w is a vector. In Fig. 2b, its element j u>x describes 
the rotational motion of the spin in the y direction and the angular velocity u in the — z 
direction. (3) From the linear spin current density j«(r, t), one can calculate (or say how 
much) the linear spin current I s flowing through a surface S (see Fig. 2d): If = JJ s dS • j s . 
However, the behavior for the angular spin current is different. From the density j w (r,t), it 
is meaningless to determine how much the angular spin current flowing through a surface 
S, because the angular spin current describes the rotational motion not the movement. On 
the other hand, one can calculate the total angular spin current 1% in a volume V from 
j u : 1% = Jff v j UJ (r,t)dV. (4) It is easy to prove that the spin currents in the present 
definitions of Eqs.(12,13) are invariant under a space coordinate transformation as well the 
gauge transformation. (5) If the system is in a steady state, j s and j u are independent of 
the time t, and 4zs(r, t) = 0. Then the spin continuity equation (10) reduces to: V • j s = j u 
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or f s dS • j s = JJJyjuidV. This means that the total linear spin current flowing out of a 
closed surface is equal to the total angular spin current enclosed. If to further consider a 
quasi one dimensional (ID) system (see Fig.2d), then one has lf — I^ = . (6) The linear 
spin current density j s = Re^^v^f?} gives both the spin direction and the direction of spin 
movement, so it completely describes the translational motion. However, the angular spin 
current density, j u = Re{^^^f} = Re{^^uj x sM/} involves the vector product of to x s, 
not the tensor us. Is it correct or sufficient to describe the rotational motion? For example, 
the rotational motion of Fig.2b with the spin s in the y direction and the angular velocity 
cD in the —z direction is different from the one in Fig. 2c in which s is in the z direction 
and lu is in the y direction, but their angular spin currents are completely the same. Shall 
we distinguish them? It turns out that the physical results produced by the above two 
rotational motions (Fig. 2b and 2c) are indeed the same. For instance, the induced electric 
field by them is identical since a spin s has only the direction but no size (see detail discussion 
below). Thus, the vector j u is sufficient to describe the rotational motion, and no tensor is 
necessary. 

Now we give an example of applying the above formulas, Eqs.(12,13), to calculate the spin 
currents. Let us consider a quasi ID quantum wire having the Rashba spin orbit coupling, 
and its Hamiltonian is: 

H = hi + v{y ' z) + S [a(x)P2 + Pxa( ^ x)] + "2^f ' (14) 

where kn(x) = a(x)m/h 2 . a(x) = for x < and x > L, and a(x) while < x < L. 
The other Rashba term —j^a x p z is neglected because z-direction is quantized^. Let ^ be a 
stationary wave function 

* (r) = V *. J tp(y,z), (15) 

2 \ e l -lo k R( x ) dx 

where (p(y,z) is the bound state wave function in the confined y and z directions, ^(r) 
represents the spin motion as shown in Fig. 3a, in which the spin moves along the x axis, 
as well the spin precession in the x-y plane in the region < x < Using Eqs.(12,13), 
the spin current densities of the wave function \l/(r) are easily obtained. There are only two 
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non-zero elements of j s (r): 

h 2 k 



j S xx(r) = '-^\ip(y,z)\ 2 cos2(f)(x), (16) 

h 2 k 

3sxy(r) = —\(p(y,z)\ 2 sm2(f)(x). (17) 



The non-zero elements of j^v) are: 



j M (r) = -^^|^)| 2 sin2<K*), (18) 



fj 2 kk 

J«»(r) = -^Ml/,*)| 8 cos20(aO, (19) 



where 0(x) = J Q X kn(x)dx. Those spin current densities confirm with the intuitive picture 
of an electron motion, precession in the x-y plane in < x < L and movement in the x 
direction (see Fig. 3a). 

In particular, in the region of x < and x > L, a(x) = kn(x) = 0, and a is a good 
quantum number, hence, j u = 0. In this case, the definition of the spin current j s is 
unambiguous. However, the spin currents are different in x < and x > L except for 
<f>(L) = rm (n = 0, ±1,±2 ...). This is clearly seen from Fig.3a. Therefore, through this 
example, one can conclude that it is sometimes impossible to define a conserved spin current. 
The example of Fig. 3a indeed exists and has been studied before^- 

Above discussion shows that the linear spin current j s = vs and the angular spin current 
ju = uj x s have clear physical meanings, representing the translational motion and the 
rotational motion (precession) respectively. They completely describe the flow of a quantum 
spin. Any physical effects of the spin currents, such as the induced electric field, can be 
expressed by j s and j u 
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IV. SPIN CURRENTS INDUCED ELECTRIC FIELDS 

Recently, theoretic studies have suggested that the (linear) spin current can induce an 
electric field _ g 21 i 22 i 23 . Can the angular spin current also induce an electric field? If so this 
gives a way of detecting the angular spin current. Following, we study this question by 
using the method of equivalent magnetic charge^. Let us consider a steady-state angular 
spin current element j u dV at the origin. Associated with the spin s, there is a magnetic 
moment (MM) m = g[J>B& — ^^-s* where //g is the Bohr magneton. Thus, corresponding 
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to juj, there is also a angular MM current j mu jdV = ^^-j^dV. From above discussions, we 
already know that j muJ (or jj) comes from the rotational motion of a MM rh (or s) (see Fig. 2b 
and 2c), and j muj = to x rh (or = uj x s). Under the method of equivalent magnetic charge, 
the MM rh is equivalent to two magnetic charges: one with magnetic charge +q located at 
Sn m and the other with — q at —5h m (see Fig.2e). n m is the unit vector of m and 5 is a tiny 
length. The angular MM current j muJ is equivalent to two magnetic charge currents: one is 
j +q = hjq8\uj\ sin 9 at the location Sn m , the other is j_ g = hjq5\uj\ sin 6* at —Sh m (see Fig.2e), 
with rij being the unit vector of j muj and 9 the angle between uj and fh. In our previous 
work,— we have given the formulae of the electric field induced by a magnetic charge current. 
The electric field induced by jm^dV can be calculated by adding the contributions from the 
two magnetic charge currents. Let 5 — > 0, and note that 2q5 —>■ \m\ and |a;||m| sin 9 = \j r 



Imu) | j 



we obtain the electric field E w generated by an element of the angular spin current j w dV: 

^^/^^/j^xl (20) 
4vr J r 6 h J r 6 

We also rewrite the electric field E s generated by an element of the linear spin current using 
the tensor j s :— 

E, = ^fvx/j,«.I (21) 

— * 

Below we emphasize three points: (i) In the large r case, the electric field E u decays as 
1/r 2 . Note that the field from a linear spin current E s goes as 1/r 3 . In fact, in terms of 
generating an electric field, the angular spin current is as effective as a magnetic charge 
current, (ii) In the steady-state case, the total electric field Et = E w + E s contains the 
property: § c Ex • dl = 0, where C is an arbitrary close contour not passing through the 
region of spin current. However, for each E w or E s , § c E w • dl or § c E s • dl can be non-zero, 
(iii) As mentioned above, a angular spin current j w may consist of different uo and s (see 
Fig. 2b and 2c). However, the resulting electric field only depends on j u = to x s. This is 
because a spin vector contains only a direction and a magnitude, but not a spatial size (i.e. 
the distance 5 approaches to zero). In the limit 5 — > 0, both magnetic charge currents j± q 
reduce to uj x fh/2 at the origin. Therefore, the overall effect of the rotational motion is 
only related to uj x m, not separately on uj and rh. Hence it is enough to describe the spin 
rotational motion by using a vector uj x s , instead of a tensor ujs. 

Due to the fact that the direction of s can change during the particle motion, the linear 
spin current density j s is not a conserved quantity. It is always interesting to uncover a 



conserved physical quantity from both theoretical and experimental points of view. Let us 
apply V» acting on two sides of Eq.(lO), we have: 

A(V.s) + V. (v'«j.-J w ) =0, (22) 

where V • j s means that V acts on the second index of j s , i.e. (V • j s )i = Ylj with 
i, j G (x, y, z). Define jv.s = V' • j s — j w , the above equation reduces into^ 

^(V«s) + V«Jv. s -=0. (23) 

This means that the current jv»s of the spin divergence is a conserved quantity in the steady 
state case. In fact, —V • s(r, t) represents an equivalent magnetic charge, so jv.s can also 
be named the magnetic charge current density. Moreover, the total electric field produced 
by '] s and j w can be rewritten as: 



E T = E s + K = ^j (V.j s -l)dVx 



r 



' x L. (24 ) 



h J r 3 

So the total electric field Et only depends on the current jv»s of the spin divergence. Note 
that Et can be measured experimentally in principle. Through the measurement of Et{t), 
js/mg can be uniquely obtained. 

In the following, let us calculate the induced electric fields at the location r = (x, y, z) 
by the spin currents in the example of Fig. 3a. Substituting the spin currents of Eqs. (16-19) 
into Eqs. (20,21) and assuming the transverse sizes of the ID wire are much smaller than 



a/?/ 2 + z 2 , the induced fields E u and E s can be obtained straightforwardly. Then the total 
electric field Et = E u + E s is: 

E T = a— V / 77 \ i^nz dx' 



m J [{x — x') 2 + y 2 + z 2 ] 3 / 2 



= aV J (V x S) . y^dx' (25) 

where V = (Hk/m, 0, 0), S = (cos 20(x'), sin 20(x'), 0), r' = (x',0,0), the constant a = 
A i o5'A i BPs/47r, and p s is the linear density of moving electrons under the bias of an external 
voltage. The total electric field Et represents the one generated by a ID wire of electric 
dipole moment p e = (0, 0, csin20(x)) at the x axis (see Fig. 3b), where c is a constant. It is 
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obvious that V x E T = 0, i.e. § c E T • dl = 0. However, in general § c • dl and § c E s • dl 

are separately non-zero. 

Finally, we estimate the magnitude of Et- We use parameters consistent with realistic 
experimental samples. Take the Rashba parameter a = 3 x 10~ n eVm (corresponding to 
/cR=l/100nm for m = 0.036m e ), p s = 10 6 /m (i.e. one moving electron per 1000 nm in 
length), and k = kp = 10 8 /m. The electric potential difference between the two points 
A and B (see Fig. 3b) is about O.OlfiV, where the positions of A and B are ^^(f; 0,0.01) 
and 2£-(f , 0, — 0.01). This value of the potential is measurable with today's technology.— 
Furthermore, with the above parameters the electric field Et at A or B is about 5V/m which 
is rather large. 

V. CONCLUSION 

In summary, we find that in order to completely describe the spin flow (including both 
classic and quantum flows), apart from the conventional spin current (or linear spin current), 
one has to introduce another quantity, the angular spin current. The angular spin current 
describes the rotational motion of the spin, and it plays a parallel role in comparison with 
the conventional linear spin current for the spin translational motion. Moreover we point 
out that the angular spin current can also induce an electric field and its E^ field scales as 
1/r 2 at large r. In addition, a conserved quantity, the current jv»s of the spin divergence, 
is discovered, and the total electric field only depends on jy.g. 
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FIG. 1: (Color online) (a) and (b) are the schematic diagram for the translational motion and the 
rotational motion of the classic vector m, respectively, (c) Schematic diagram for a classic vector 
flow. 

FIG. 2: (Color online) (a) The linear spin current element j s ,xy (b) and (c) The angular spin 
current element ju jX . (d) The spin current in a quasi ID quantum wire, (e) The currents of two 
magnetic charges that are equivalent to a angular MM current. 

FIG. 3: (Color online) (a) Schematic diagram for the spin moving along the x axis, with the spin 
precession (rotational motion) in the x-y plane while < x < L. (b) A ID wire of electric dipole 
moment p e . This configuration will generate an electric field equivalent to the field from the spin 
currents in (a). 
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